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ABSTRACT 

We present the results of a series of gas dynamical cosmological simulations of the 
formation of individual massive field galaxies in the standard concordance ACDM and 
in a AWDM cosmology each with f2 =0.3 and A o =0.7. Two high resolution simula- 
tions (2 x 50 3 gas and dark matter particles) have been performed and investigated 
in detail. The gas component was represented by Smooth Particle Hydrodynamics 
(SPH) and a simple star formation algorithm was applied. The galaxies form in an 
initial burst of star formation followed by accretion of small satellites. They do not 
experience a major merger. The simulated galaxies are old (« 10 Gyrs) hot stellar 
systems with masses of w 1.7 x 1O 11 M0. Baryonic matter dominates the mass in the 
luminous part of the galaxies up to « 5 effective radii. The projected properties of 
the galaxies have been investigated in detail: The ACDM galaxy is a slowly rotating 
((w/er) max = 0.2) spheroidal stellar system (E2) with predominantly disky isophotes. 
The line-of-sight velocity distributions (LOSVDs) deviate from Gaussian shape and /13 
is anticorrelated with vi os . The corresponding AWDM galaxy is more elongated (E3 - 
E4) and rotates faster ((v/er) max = 0.6). The anisotropy parameter (v/a)* is close to 
unity indicating isotropic velocity dispersions. There is no clear indication for isopho- 
tal deviations from elliptical shape and the projected LOSVDs do not show correlated 
higher order deviations from Gaussian shape. Within the uncertainties of M/L both 
galaxies follow the Fundamental Plane. We conclude that the properties of the two 
galaxies simulated in the ACDM and AWDM cosmology are in good agreement with 
observations of intermediate mass elliptical or SO galaxies. Our conclusion differs from 
Meza et al. (2003), who find, from a similar simulation, a much more concentrated 
galaxy than is generally observed. The differences in our findings may either be the 
result of differences in the star formation algorithms or due to the different merger 
history of the galaxies. 
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1 INTRODUCTION 

The concordance ACDM paradigm (a cold dark matter cos- 
mology with the addition of a cosmological constant) ap- 
pears to provide an excellent fit to astronomical observations 
on scales large compared to the sizes of individual galaxies 
(Percival et al. 2002; Spergel et al. 2003). However there 
are some indications that the standard model may have too 
much power on small scales to be consistent with observa- 
tions. For example, Klypin et al. (1999) show that numeri- 
cal simulations predict a larger number of Galactic satellites 
than observed, though Benson et al. (2002) argue that this 
problem can be solved by the suppression of dwarf galaxy 
formation in a photoionized inter-galactic medium. A sec- 



ond problem relates to the steep cusps found in the centres 
of simulated dark matter halos (see Moore et al. 1999; Eke 
et al. 2001; Binney & Evans 2001), which appear to be incon- 
sistent with the dark matter distributions inferred in dwarf 
galaxies (see e.g. Dutton et al. 2003). It is not yet clear 
whether this discrepancy requires a revision of the ACDM 
model. For example, Ricotti (2003) argues that low mass 
haloes in the CDM model may have less cuspy profiles than 
higher mass haloes, though this result is disputed by Colin 
et al. (2003). 

It is also not yet clear whether the properties of real 
galaxies can be explained by the ACDM model. Does the 
concordance model produce galaxies of the right masses and 
sizes at the right epochs? In fact there are some indica- 
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tions from galaxy formation that the concordance model 
may have too much small scale power. For example, it has 
proved difficult to make realistic disk galaxies in numeri- 
cal simulations of the CDM model incorporating gas dy- 
namics. In most simulations, the disk systems that form are 
smaller, denser and have much lower angular momenta than 
real disk systems (see Navarro & White 1994, Navarro et al. 
1995; Navarro & Steinmetz 1997, Steinmetz & Navarro 1999; 
Wright 2002). More acceptable fits to real disk systems can 
be found if heuristic prescriptions modelling stellar feedback 
are included in the simulations (Weil et al. 1998, Sommer- 
Larsen et al. 2002,Governato et al. 2002, Abadi et al. 2003). 
However, even in these simulations, the disk systems typi- 
cally contain denser and more massive bulges than the vast 
majority of real disk galaxies. 

Most of the previous work on the formation of individ- 
ual galaxies from cosmological initial conditions has focused 
on the formation of disk galaxies. The formation of indi- 
vidual elliptical galaxies has not been investigated as exten- 
sively. This seems surprising as giant elliptical galaxies are 
the oldest and most massive stellar systems in the Universe 
and probably contribute over 50% of the total stellar mass 
if we include the stars in the bulges of SO, Sa and Sb spi- 
rals. Although their internal kinematics can be very complex 
the major component of the stellar population in ellipticals 
is old and homogeneous. They are therefore good probes of 
galaxy assembly, star formation and metal enrichment in the 
early universe (see e.g. Thomas et al. 2002). Furthermore, 
the giant ellipticals follow simple scaling relations, the Fun- 
damental Plane being the most important (see e.g. Bender 
et al. 1992). These simple scaling relations should arise nat- 
urally from the correct cosmological model. 

Despite their complex kinematics, it has become evi- 
dent over the last 15 years that observed giant ellipticals 
show detailed photometric and kinematic properties that 
correlate with their luminosity. Massive giant ellipticals are 
slowly rotating, flattened by anisotropic velocity dispersions 
and show predominantly boxy isophotes. Lower mass giant 
ellipticals have disky isophotes and are flattened by rota- 
tion (Bender et al. 1988). These low mass ellipticals most 
likely contain weak disk components (Rix & White 1990). 
The fact that boxy ellipticals, in contrast to disky ellipti- 
cals, show strong radio and X-ray emission (Bender et al. 
1989) and have flat density cores (Faber et al. 1997) might 
indicate that they formed either by a different process or in 
a different environment. 

How and when giant ellipticals have formed is still 
poorly understood. According to the "merger hypothesis" 
early type galaxies formed by mergers of disk galaxies. Ide- 
alised models of binary mergers of disk galaxies (with and 
without gas and star formation) and multiple mergers have 
been investigated in great detail by several authors (e.g. 
Negroponte & White 1983; Hernquist 1992; Barnes 1988; 
Barnes & Hernquist 1996; Mihos & Hernquist 1996; Weil 
& Hernquist 1996; Springel 2000). Those simulations - the 
recent ones with high numerical resolution - are useful in 
understanding detailed internal processes e.g. gas inflow to 
the centre (Mihos & Hernquist 1996). They are also capable 
of explaining the origin of fine structure in individual ellipti- 
cals. For example, a large study of collisionless disk mergers 
by Naab & Burkert (2003) showed that binary disk mergers 
can successfully reproduce global kinematic and photomet- 



ric properties of low and intermediate mass giant ellipticals. 
The formation of faint embedded disks that are observed in 
these galaxies can be explained if gas was present in progen- 
itor galaxies (Naab & Burkert 2001a, b; Barnes 2002). 

Despite these successes, binary merger simulations suf- 
fer from certain limitations. In particular they use approx- 
imate equilibrium models of present day spiral galaxies 
as progenitors rather than self-consistently calculating the 
properties of the progenitors "ab initio" from realistic cos- 
mological initial conditions. This is a serious limitation, since 
it is unlikely that the high redshift progenitors of ellipticals 
really resembled present day spirals. Added to the model 
uncertainties there are also a large number of degrees of 
freedom in the initial conditions, e.g. the geometries of the 
orbits, halo profiles, bulge masses, bulge rotation, gas con- 
tent, gas distribution, disk sizes etc. are all adjustable pa- 
rameters. Although there have been attempts to survey e.g. 
different halo profiles (Dubinski et al. 1996) or disk spin ori- 
entations (Naab & Burkert 2003), these parameter surveys 
are evidently incomplete. In addition, questions regarding a 
self consistent evolution of stellar populations are extremely 
difficult to address. Specifically current generation ellipti- 
cals are far too red, metal rich and old to have formed via 
mergers of systems similar to current epoch spirals. Merger 
simulations have therefore failed, so far, to explain the origin 
of global scaling relations like the color-magnitude relation 
or, more generally, the fundamental plane. They have, how- 
ever, proved very useful in developing an understanding of 
the detailed internal merger dynamics. 

The best way to overcome these problems is via high 
resolution simulations of individual elliptical galaxies from 
realistic cosmological initial conditions. The initial condi- 
tions are constrained by the cosmological model alone and 
the subsequent evolution is governed solely by the numerical 
resolution and accuracy of the physics that is implemented 
in the simulation. Once a sufficient number of individual el- 
lipticals over the whole mass spectrum have been simulated, 
it should be possible (if the cosmological model is correct) 
to explain the origin of the global scaling relations and the 
detailed properties of individual galaxies at different lumi- 
nosities. 

A first attempt in this direction has been made by Meza 
et al. 2003. These authors followed the formation of a sin- 
gle spheroidal galaxy in a ACDM cosmology. The kinematic 
properties of their galaxy resembled a rotationally supported 
giant elliptical. However, the effective radius of the simulated 
galaxy was a factor of 10 smaller than for observed ellipti- 
cals at the same brightness. This galaxy is therefore much 
too compact to be consistent with observations. 

Does this result imply a problem in forming elliptical 
systems in the ACDM model? Two lines of investigation 
are suggested. First, the cosmological model adopted may 
be correct but the physical treatment may be inaccurate. 
Specifically, feedback from some early star formation into 
the shallow potential wells in the small halos existent at 
those times, may so efficiently blow out other gas as to re- 
duce early star formation effectively (Dekel & Silk 1986, Ef- 
stathiou 2000, Springel & Hernquist 2003 and Nagamine, 
Cen and Ostriker (in preparation) are also exploring this 
possibility). The inclusion of stellar feedback would reduce 
the number of low mass galaxies, but not the number of 
low mass halos. It would also significantly reduce the stellar 
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density in the centres of systems, but would not reduce the 
dark matter density by very much 1 . 

If, however, it is not obvious that the discrepancy found 
by Meza et al. 2003 can be cured by a better treatment of 
the physics. If the evidence for low dark matter densities in 
the centres of galaxies is taken seriously (Binney & Evans 
2001), then we may want to consider a more radical solution. 

An alternative is to reduce the small scale power in 
the dark matter density fluctuations. This can be achieved 
in various ways. First, the spectral index, n, could be suf- 
ficiently small so that after normalization to the WMAP 
amplitude, and extrapolation to the small wavelengths rel- 
evant to galaxy formation, the amplitude is low enough to 
significantly reduce early star formation. The WMAP anal- 
ysis (Spergel et al. 2003) combined with 2dFRGS and su- 
pernova data in fact indicated that n may be as small as 
0.93 ± 0.03, and this may alleviate some of the purported 
difficulties of the concordance ACDM model. We will return 
to this possibility in later work. However, there are signifi- 
cant limits on the value of the spectral index n, since infor- 
mation on the cluster length scale (~ 8/i _1 Mpcj, which is 
intermediate between the WMAP scale and the galaxy for- 
mation scale seems to require a relatively high normalization 
(Bahcall et al. 2002). This would not permit a constant n 
solution with n much less than 0.95. Furthermore, the high 
optical depth for electron scattering becomes very difficult 
to achieve with a low spectral index (Cen 2003; Chiu et al. 
2003). 

Another possibility is to achieve low power on small 
scales by some form of cut-off in the power spectrum. An 
example of this is the Warm Dark Matter model, where a 
finite but quite small initial "thermal" velocity dispersion 
sharply truncates the power as (k/k cu t)~ 10 above some wave 
number scale, k cu t- Bode, Ostriker & Turok (2001) show that 



of a high optical depth to the surface of last scattering, r es , 
rule out WDM? The answer is ambiguous. An examination 
of Fig. 5 of Spergel et al. (2003) indicates for n = 0.95, 
r es may be 0.05 at la level and as small as 0.01 at the 
2<t level. Further would is needed (but see Yoshida et al. 
2003 for a counter argument) to set a limit on m x based 
on WMAP. Detailed work in progress by Ricotti & Ostriker 
(2003) indicates that, if other parameters are held constant, 
a WDM model with m x = 1.25 KeV leads to only a ~ 10% 
reduction in r es . 

In this paper we aim to investigate formation and evo- 
lution of intermediate mass giant galaxies in the ACDM 
model, and to quantify the effect of reducing the power 
at small scales by studying a AWDM model. The paper is 
organised as follows: Section 2, summarizes the simulation 
code and describes how ACDM and AWDM initial condi- 
tions were generated . The results of two high resolution 
simulations in the AWDM and ACDM universe and a com- 
parison of their global properties with a set of low resolution 
simulations are described in Section 3. In Section 4 we com- 
pare the internal properties of the two high resolution simu- 
lations in detail with observations of giant elliptical galaxies. 
Section 5 contains a summary and our conclusions. 



2 INITIAL CONDITIONS AND SIMULATIONS 

We performed two sets of cosmological simulations using 
ACDM and AWDM initial conditions: 5 simulations with 
2 x 34 3 gas and dark matter particles and 1 simulation with 
2 x 50 3 gas and dark matter particles. The lower resolu- 
tion simulations were used to investigate the differences in 
the global properties of the stellar systems formed in each 
cosmology. The higher resolution simulations were used to 
investigate the internal properties of the simulated galaxies. 
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where the warm dark matter particle has mass, m x and the 
density in WDM is represented by Q x , at z = 0. The addi- 
tional velocity dispersions are distributed as 



f(v) = (e^ + iy 1 
with 



(2a) 



Barkana, Haiman & Ostriker (2001) used the abundance of 
small column density lines in the Lyman alpha forest to 
limit m x > 0.75 keV . Other work suggests similar limits of 
m x > 1 keV , hence, in this paper, we will investigate WDM 
with m x — 1 keV . 

One may ask at this point if the WMAP observations 



A small reduction would occur since concentration of baryonic 
material forces a moderate increase of the central dark matter 
density by purely gravitational processes over what would have 
been the case without efficient baryonic cooling. 



2.1 Power Spectra & Cosmological Parameters 

The initial conditions of the ACDM and AWDM simulations 
assumed scale-invariant adiabatic fluctuations. The post- 
recombination power spectrum for both cosmologies is based 
on the parameterisation of Efstathiou, Bond & White (1992) 
with T=0.2 and the transfer function 



T(k) = 



where 



{1 + [ak + (6fc) 3 /2 + i^yYYI" ' 



•-(¥)■'-(¥)•»(¥) 

with units of h 1 Mpc and 
v = 1.13 and V = 0.2. 



(3a) 



(3b) 



(3c) 



The power spectrum of the dark matter fluctuations in 
the ACDM cosmology is given by 



P(k) oc T 2 (k)k. 



(4) 



The WDM spectrum is truncated at small scales as de- 
scribed by Bode, Ostriker & Turok (2001) through an addi- 
tional factor in the transfer function of the CDM spectrum 
such that 
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Figure 1. (a) The particle distribution of the AWDM simulation initial conditions, (b) A comparison between the input 3-D power 
spectrum for the ACDM cosmology and a warm dark matter power spectrum with particle mass, m x = 1.0 keV. (c) The calculated 
3-D power spectrum from the initial conditions of the ACDM and AWDM simulations. The dotted line is the Nyquist frequency of the 
particle grid in the case of the 50 3 particle simulations and is also the wavenumber at which noise dominates the spectrum as a result 
of the discrete nature of the initial conditions. We note that the AWDM power spectrum has a cut off which occurs at wave-numbers 
lower than this numerical limit. 



TwDM(fc) = [i+&)r (5) 

with k cu t given by Equation 1. Here, m x is the mass of 
the warm dark matter particle which is set to 1 keV as 
described in the Introduction. The AWDM power spectrum 
P(k) therefore takes the following form: 

P(k) « TvV DM (fc)T 2 (fc)fc. (6) 

The amplitude of the mass fluctuations in both cos- 
mologies is normalised so that the rms mass fluctuation in 
spheres of radius 8fe _1 Mpc is <T8=0.86. This normalization 
was chosen to match the present day abundances of rich 
galaxy clusters in this cosmology (Eke et al. 1996) and is 
somewhat higher than the value deduced from more recent 
analyses (Pierpaoli et al. 2002, Seljak 2001 and Viana et al. 
2002). However, the value of as = 0.86 is consistent with the 
WMAP determination and with weak lensing constraints 
(Refrigier 2003), for our adopted value of the spectral index 
n = 1. 

The baryonic fraction is set to ft = f2(,/fi m =0.17, con- 
sistent with the predictions of primordial nucleosynthesis 
(Buries & Tytler 1998) and WMAP (f b = 0.171 ±0.025) for 
a Hubble parameter of h=0.65 2 .. Figure 1(a) shows a 2-d 

2 h is defined such that _ffo=100h kms _1 Mpc _1 . 



projection of a the initial particle positions for our 50 3 sim- 
ulation in the AWDM cosmology. Figs. 1(b) and (c) show 
the power spectra of the ACDM and AWDM initial condi- 
tions calculated from the particle distributions at z = 24. 
The dotted line in Figure 1 (c) shows the limit of our spa- 
tial resolution in the 50 3 simulations. At scales smaller than 
this limit (to the right of the dotted line) there is no input 
power and the Poisson noise dominates. 

Velocity dispersions were added to the initial conditions 
of the warm dark matter model according to equations (2a 
& 2b). Additional velocity dispersions smear out the small- 
scale perturbations through free-streaming. We have per- 
formed simulations both with and without additional veloc- 
ity dispersions and find only small differences at the very 
centres of the dark matter haloes (slightly less steep cusps 
when the velocity dispersions were included). The effects of 
initial velocity dispersions are unimportant for the simula- 
tions described in this paper. 

We used identical values of S7q=0.3, Ao=0.7 for both 
the AWDM and the ACDM simulations. The mass of the 
warm dark matter corresponds to a mass cut-off of 

/ O \ -°- 45 / h \ 3 9 / m \ ~ 3 - 45 

•-'"^(n) (4) im) < 7 » 

in the AWDM simulations. The mass resolution of the 50 3 
AWDM simulation is already sufficient to correctly repre- 
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Simulation 


Number 


e g (kpc) 


«d (kpc) 


v c (kms 1 ) 


m dm (!0 6 M e ) 


m g (1O 6 M ) 


ACDM 34 3 


5 


0.8 


1.4 


126 - 167 


116 - 208 


29 - 52 


ACDM 50 3 


1 


0.5 


0.85 


165 


66 


16 


AWDM 34 3 


5 


0.8 


1.4 


125 - 169 


116 - 208 


29 - 52 


AWDM 50 3 


1 


0.5 


0.85 


162 


66 


16 



Table 1. Summary of some of the properties of the ACDM SPH simulations that have been re-simulated in the AWDM cosmology. The 
first column gives the number of gas particles within the high resolution cube (equal to the number of dark matter particles within this 
volume). The second column lists the number of simulations at each mass resolution. The third and fourth columns list the Plummer 
softening parameter for the gas and dark matter particles respectively. The fifth column lists the range of the circular speeds of the dark 
matter haloes (determined within the virial radius at z = from the low resolution AP 3 M simulations) . The sixth column lists the range 
of gas particle masses. 



sent the input power spectrum above this mass cut-off (see 
Table 2) , thus the distribution of dark matter would not be 
expected to change in a simulation with still higher resolu- 
tion. However, this is not true of the ACDM simulations, 
which have no mass cut-off. Higher mass resolution in the 
ACDM would be expected to produce earlier collapse of 
smaller structures. 



2.2 Halo Selection 

The initial conditions for our simulations were generated 
in a two stage process. A low resolution dark matter only 
simulation of a large computational volume was run until 
the present day. Haloes were selected from the final output 
of the simulation to be run at higher resolution including 
Smooth Particle Hydrodynamics (SPH) to represent a dissi- 
pative gas component. The re-simulation of these haloes at 
higher resolution requires the addition of short wavelength 
fluctuations missing in the low resolution simulations. In ad- 
dition particles at large distances are averaged hierarchically 
to represent the tidal field of the entire computational vol- 
ume by a manageable number of particles. A full description 
of the methods used to generate our initial conditions can 
be found in Weil, Eke & Efstathiou (1998) (see also Navarro, 
Frenk & White 1997). Here we summarise the specific pa- 
rameters used for the simulations in this paper. 

We ran a low resolution dark matter simulation us- 
ing the Adaptive Particle-Particle, Particle-Mesh (AP 3 M) 
N-body code of Couchman (1991). For both cosmologies, a 
cube of Lbox = 50 Mpc containing 128 3 particles was evolved 
from a redshift of z = 24 to the present day. The force law 
used in the code was a Plummer law with a gravitational 
softening fixed in comoving coordinates at 7 kpc . Identi- 
cal random phases were used in both cosmologies so that 
the same haloes could be simulated at higher resolution en- 
abling direct comparisons to be made. We used a spherical 
over-density group-finding algorithm (Lacey & Cole 1994) to 
identify virialised haloes at z — with masses in the range 
of 2 x 10 12 M Q < M halo < 7 x 10 12 in a low density environ- 
ment such that the nearest halo with a mass greater than 
7.1 x lO n M is over 1 Mpc away. If applied to a 50 3 Mpc 3 
volume of the real universe our selection procedure would 
identify intermediate-mass giant field ellipticals or SO galax- 
ies systems similar to the Sombrero galaxy (M104) rather 
than late type spirals which are of lower mass, or very bright 
giant ellipticals which are found in high density regions. 



Each target halo was re-simulated at a higher resolution 
with the inclusion of a gas component. We increased the par- 
ticle number within a cubic volume at redshift z = 24 con- 
taining all particles that end up within the virialised region 
of the halo at z = 0. Additional short wavelength pertur- 
bations were included to account for the missing small-scale 
power below the Nyquist frequency of the low resolution 
simulation. The size of the high resolution cube was of the 
order of Lb ox = 5kpc. The tidal forces from particles outside 
the high resolution cube were approximated by increasingly 
massive particles in 5 nested layers of lower and lower res- 
olution. Gas particles with masses chosen to the match the 
baryonic fraction were added in the high-resolution regions, 
in the same number and at the same positions as the dark 
matter particles. Outside the high-resolution regions only 
dark matter particles were included. 

In total we ran 6 simulations simulations from z = 24 
to z = for each of the AWDM and ACDM cosmologies: 
Five simulations with 2 x 34 3 particles (labelled (a) to (e)) 
and one simulation with 2 x 50 3 particles within the high 
resolution cube . The simulations will be referred to as the 
34 3 and 50 3 simulations hereafter. At the mass and spatial 
resolution of the 34 3 runs, it is not possible to resolve the 
internal structure of a forming galaxy and they are used in- 
stead to investigate variations in their global properties. The 
50 3 simulations are used to investigate the internal proper- 
ties of the galaxies in detail. 

The simulations were performed with the GRAPESPH 
code as outlined in Weil et al. (1998) using 5 GRAPE-3A 
boards (Sugimoto et al. 1990) connected to a Sun Ultra-2 
workstation. The evolution of the gas component was fol- 
lowed with SPH including radiative cooling and a simple 
star formation algorithm: Each gas particle that remains in 
a collapsing region with a density p cr ; t >7xl0~ 23 kgm -3 (see 
Navarro & White 1993; Weil et al. 1998) for longer than a lo- 
cal dynamical time was converted into a star particle. Stellar 
feedback has been included in some simulations, using some- 
what 'ad hoc' rules (Sommer-Larsen et al. 2002; Meza et al. 
2003; Springel & Hernquist 2003). However, it is clear that 
these representations are far from realistic. We have there- 
fore chosen to make our simulations as simple as possible 
and hence stellar feedback is ignored. This makes it easier to 
compare simulations done at different mass resolutions, and 
to understand if the properties of the final stellar systems 
depend on whether the dark matter is warm or cold. The- 
oretical arguments (see e.g. Dekel & Silk 1986; Efstathiou 
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a 




b 




c 




d 




e 




ACDM 


AWDM 


ACDM 


AWDM 


ACDM 


AWDM 


ACDM 


AWDM 


ACDM 


AWDM 


Rmr/kpc 


350 


360 


330 


330 


280 


280 


280 


280 


300 


310 


M O Af(R^r)/lO n M 


21 


22 


17 


17 


10 


10 


10 


11 


13 


14 


loglO (jh (Rvir)) 


3.1 


3.1 


3.2 


3.1 


3.3 


3.4 


3.3 


3.3 


3.2 


3.3 


M 9a! /1O 1O M 


8.0 


8.9 


7.3 


7.9 


7.4 


6.1 


5.5 


6.9 


6.3 


7.2 


!ogio (i*/kpc kms -1 ) 


1.7 


2.3 


2.1 


1.2 


1.4 


1.5 


2.6 


2.7 


2.0 


2.2 



Table 2. The virial radii, virial masses, galaxy masses and specific angular momenta of the galaxies (measured in kms -1 kpc ) and dark 
matter haloes for each of our five, 34 3 simulations in the ACDM and AWDM simulations at z=0. The difference between the cold and 
warm dark matter simulations are not statistically significant at the gross, integrated level shown in this table. 
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Figure 2. The number of groups (averaged over the 5 simula- 
tions) of a given mass in the 34 3 AWDM (solid) and ACDM 
(dashed) simulations for both dark matter (thick) and stars (thin) 
at z = 0. 



2000) suggest that stellar feedback is less important in high 
velocity dispersion systems than in low velocity dispersion 
systems, though since galaxies form hierarchically, feedback 
processes may have a significant effect on the evolution of 
high velocity dispersion systems at early times. 



3 ACDM VS. AWDM 



3.1 34 3 Simulations 

For the AWDM cosmology investigated in this paper the 
formation of halos with masses M haio < M cut — 1 x 1O 1O M0 
should be suppressed because the initial power spectrum is 
truncated below this mass-scale. To quantify this effect, we 
compared the mass spectrum of satellite halos in the two 
cosmologies at the present day. A friends-of-friends group 
finding algorithm with a linking length of 0.8 kpc was used 
to identify dark matter haloes and their stellar systems for 
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Lookback time (Gyrs) 



Figure 3. Normalised stellar mass of the central object versus 
lookback time, averaged over the 5 simulations in the ACDM 
(dashed) and AWDM simulations (solid). The error bars indicate 
the one sigma error. On average the ACDM galaxy assembles its 
stars earlier than the AWDM galaxy. 



the 34 3 simulations at z = . The massive central galaxy 
was excluded. 

Figure 2 shows the mass distributions of the dark 
matter and stellar systems averaged over the 5 simula- 
tions in each cosmology As expected, the ACDM cosmol- 
ogy produces many more dark matter haloes with masses 
< 1 x 1O 1O M0. There is, however, only a marginal difference 
in the mass distribution of the stellar systems. The central 
galaxies in the 34 3 simulations have halo masses in the range 
of 1.0 x 1O 12 M to 2.2 x 1O 12 M , well above the cut-off mass 
of the AWDM cosmology. Table 2 summarizes some of the 
global properties of the ACDM galaxies are their AWDM 
counterpart. The virial radii, masses and the specific angu- 
lar momenta of the galaxies and their halos do not reveal 
any large differences between the cold and warm dark mat- 
ter simulations. 

A comparison of the assembly history of the stellar mass 
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Figure 4. 

at z=0. 



x (kpc) 

Three orthogonal projections of the star particle positions of the (a) ACDM and (b) AWDM galaxy in the 50 3 simulations, 



for the central ACDM and AWDM galaxies averaged over 
the five 34 3 simulations is shown in Figure 3. The mean 
values and standard deviations are shown. The redshifts at 
which 1/10 and 1/2 and 9/10 of the present day mass in 
stars has been assembled are listed in Table 3. Again, there 
is no obvious difference between the ACDM and the AWDM 
cosmologies. At z = the stellar systems in both cosmologies 
have mean ages of about 10 Gyrs. The masses and ages 
of the galaxies are thus comparable to those of early type 
galaxies. In conclusion, despite differences in the formation 
of low mass dark matter haloes, the massive central galaxies 
formed in the two cosmologies resemble early type galaxies 
and have indistinguishable formation histories and global 
properties. The resolution of the 34 3 simulations, however, 
is too low to investigate the internal properties of the main 
stellar systems in detail. We therefore now turn to the higher 
resolution 50 3 simulations. 
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Figure 5. The radial density profile for the galaxies in our 50 3 
(a) AWDM simulation and (b) ACDM simulation. Two exponen- 
tial fits are also plotted in the figure to represent an inner and 
outer component. The combined fit is also shown. The best fitting 
exponentials have been calculated according to Eqn. 8. 



3.2 50 3 simulations 

We have re-simulated halo 'a' (see Table 2) using 2 x 50 3 
particles within the high resolution region. Figure 4 shows 
the three orthogonal projections along the principal axes of 
the moment of inertia tensor of the central stellar systems in 
the ACDM and AWDM cosmologies at z = 0. The ACDM 
simulation produces a spheroidal stellar distribution. The 
AWDM galaxy on the other hand has a slightly more flat- 
tened morphology. We assume for a moment that the galax- 
ies are two component systems and perform a decomposition 



of the surface mass density profile into an inner and outer 
component of the form: 



lo Sio 



H{r) 



/to 



= logio 



a exp 



(8) 



where the constants a, b, c and d were determined using a 
non-linear fitting algorithm. The surface density profiles and 
the fits for the two galaxies are shown in Figure 5. The stel- 
lar system that forms in the AWDM simulations has a more 
extended outer component. More interestingly, however, we 
see that the inner component of the galaxy in the AWDM 
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simulation is more extended than in the ACDM galaxy. To 
test for evidence of whether the inner and outer components 
have kinematic counterparts, we performed a kinematic de- 
composition of the stars. We used a weak criterion to sepa- 
rate a kinematically cold component from a hot component. 
A star was tagged as a member of a cold component if 



> 0.7« c 



ve,i > 



(Vrand) 



9/ 2 vl/2 

+ (v$,bin 



ve, 



(9a) 
(9b) 

(9c) 



where v c is the average circular velocity of the radial bin 
containing the star, v r ,i is the radial velocity of the star, 
v$,i is the rotational velocity of the star, v z ,i is the velocity 
along the minor axis defined by the moment of inertia tensor 
and v g>bin is the average rotational velocity of the radial bin 
containing the star. Both galaxies are essentially hot stellar 
systems with no significant indication of a dynamically cold 
sub- component. The AWDM galaxy contains a small num- 
ber of particles which obey the criteria given in equations 
9a - 9c which are, however, distributed in a thick disk-like 
structure. 

To quantify whether the flattening of the AWDM 
galaxy is caused by rotation, we have calculated the dimen- 
sionless spin parameter, A, for the dark matter halos and the 
stellar systems where A is defined as 
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Figure 6. The radial profile of the stellar specific angular mo- 
mentum as a fraction of the specific angular momentum calcu- 
lated from the circular velocity. The main stellar system in the 
ACDM (dashed) cosmology has only weak rotational support. 
The AWDM (solid) galaxy has more angular momentum in the 
outer regions. 



A = 



J\E\ 



1/2 



GM 5 / 2 • < 10 > 

and J is the angular momentum, E is the total energy and 
M is the total mass of the system (see e.g. Fall & Efstathiou 
(1980)). To compute (10) only bound halo particles inside 
the virial radius and bound stars within the outer radius of 
the luminous part of the galaxy are included (defined as the 
radius at which the stellar density of the galaxy averaged in 
spherical shells becomes equal to the background level). The 
stellar system is treated as if it were virialised and so \E\ in 
Eqn. 10 is replaced by the kinetic energy of the stars |T*|). 
In the ACDM cosmology the halo has Ah=0.041 and the 
stars have A* =0.053, in the AWDM cosmology the halo has 
a slightly lower value of A.sr=0.032 whereas the stars have 
higher value of A» =0.228. Fig. 6 shows the specific angular 
momentum profiles for the ACDM and the AWDM galaxy. 
Neither cosmology produces a purely rotationally supported 
stellar system (j*{r)/j c (r) ~ 1, where j c {r) = rv c (r) is 
the specific angular momentum computed from the circular 
speed). However, the stellar system of the AWDM galaxy 
has a higher angular momentum than the ACDM galaxy 
except in the central regions r < lkpc. 

Figure 7 shows the time evolution of the average specific 
angular momentum of the dark matter and stars in both 
cosmologies. At any given redshift the dark matter halo of 
the ACDM galaxy has more angular momentum than the 
AWDM halo, but they both show a very similar evolution 
with time. In contrast, the angular momentum of the stars 
of the ACDM galaxy does not change significantly after a 
redshift of z — 4 whereas the stars of the AWDM galaxy gain 
further angular momentum until z = 1 with little evolution 
thereafter. 

We have investigated the merger histories of the galaxies 
in more detail. A friends-of-friends algorithm with a linking 
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Figure 7. The evolution of the specific angular momentum of 
stars (thick) and dark matter (thin) for the ACDM (dashed) and 
AWDM (solid) galaxy. 



length of 2.0 kpc (> 2e) was used to identify satellite galax- 
ies consisting of more than 20 stellar particles (3.2 x 1O 8 M0) 
at each redshift. A satellite was then defined to have merged 
with the central galaxy at a redshift at which its particles 
were first identified as members of the central galaxy. Figure 
8 shows the stellar mass ratio of every identified merger as a 
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Figure 8. Ratio of the stellar mass of satellites that merge with 
the central galaxy to the stellar mass of the galaxy as a function 
of lookback time for both the ACDM (open squares) and AWDM 
(filled triangles) cosmology. Every identified satellite with a mass 
larger than 3.2 X 10 18 M Q is plotted. 



function of lookback time. The total number of mergers with 
M* a t > 10 9 M is 58 in the AWDM simulation compared 
with 96 in the ACDM simulation. The results show that at 
almost every redshift the ACDM galaxy experiences more 
individual minor mergers with higher mass ratios in stars 
than the corresponding AWDM galaxy. The mass ratio for 
an individual merger never exceeds 10:1. Low accretion rates 
like this are expected as the initial conditions were selected 
from low density environments. In addition to the frequency 
of mergers, the average internal composition of the merg- 
ing satellites changes with time and is different in the two 
cosmologies (Figure 9) . At redshifts z < 2 the ACDM satel- 
lites always have a higher star-to-gas ratio. Interestingly, the 
only redshift range where the satellites have more mass in 
gas than in stars is 2 < z < 1. This corresponds to the period 
where the galaxies reach their second peak in the gas-to-star 
mass ratios as shown in Figure 10. In this figure we plot the 
time evolution of the composition of the central galaxy. The 
mass-ratios of gas to stars (top) and dark matter to stars 
(bottom) in the AWDM and ACDM galaxies are plotted. 
The dark matter fractions vary only slightly between the 
two cosmologies. However, between a redshift of z — 4 and 
z — 1 the AWDM galaxy is significantly more gas rich than 
the ACDM galaxy. 

To investigate where the stars that end up in the final 
galaxies have formed, we have calculated the stellar mass 
accumulated by the galaxies through mergers and compared 
it to the cumulative mass of stars that formed inside the 
galaxies (Figure 11). Below a redshift of two, when both 
galaxies have already assembled ?a 50% of their final stellar 
mass (see Figure 16), about 30% of all present day stars 
in the ACDM galaxy have been accreted by mergers and 



Figure 9. The star to gas mass-ratios of merged satellites as 
a function of lookback time for the AWDM (solid) and ACDM 
(dashed) simulations. Satellites merging at z < 1 are in general 
more gas rich than ACDM satellites. 
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Figure 10. The time evolution of the mass-ratio of gas to stars 
(top) and dark matter to stars (bottom) of the AWDM (solid) and 
ACDM (dashed) galaxy within 20 kpc. Throughout the evolution 
the AWDM galaxy is more gas rich than the ACDM galaxy. 
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Figure 11. Top: The cumulative stellar mass accreted in mergers 
as a function of lookback time for the AWDM galaxy (solid) and 
the ACDM galaxy (dashed). Bottom: The evolution of the stellar 
mass formed within the galaxies. The masses include material 
inside a radius of 20 kpc. After z = 2 the galaxy in the AWDM 
simulation accretes fewer stars from mergers and interactions and 
forms the majority of its new stars within the galaxy itself. The 
ACDM galaxy, however, accretes nearly half its mass in stellar 
mergers. 



Figure 12. The average age of all the stars within the galaxy 
in the AWDM (filled squares) and ACDM (open squares) cos- 
mologies. Linear least squared fits are plotted for the average 
age of stars formed in the galaxy (thin) and the average age of 
stars accreted in mergers (thick) for ACDM (dashed) and AWDM 
(solid). The dotted line represents the age of the universe. The 
stellar populations of the mergers in the AWDM simulation are 
on average younger than the corresponding ACDM populations. 
At z < 1 there is a difference of approximately 1 Gyr in age. 



~ 20% of the stars have formed inside the galaxy. Over the 
same period, the ACDM galaxy has accreted only « 10% 
of its present day stars and « 40% have formed within the 
galaxy. 

The influence of the different accretion histories on the 
mean ages of the stellar population of the galaxies is shown 
in Figure 12. In each redshift bin we show the average age 
of the total stellar population for the ACDM (open squares) 
and the AWDM (filled squares) galaxies. In addition, the 
total stellar population was divided into stars that have 
merged with the galaxy (thick dashed and solid lines) and 
stars that have formed in the galaxy (thin dashed and solid 
lines). At all redshifts the ACDM galaxy accretes stars that 
are of a similar age to the galaxy itself. The accreted stars 
have formed in small halos at the same time as the stars 
within the main galaxy. In contrast, the stars accreted in 
the AWDM universe below z — 0.5 are on average 1 Gyr 
younger than the stars within the galaxy. In total, however, 
the number of accreted stars is so small that there is almost 
no effect on the mean age of the AWDM galaxy. 

Fig. 13 shows the star formation rate of the main galaxy 
in the AWDM and ACDM simulations. The star formation 
rate of the ACDM galaxy can well be approximated by an 
exponential, sfr oc exp (— (t — to)/r), with to = 13.5 Gyrs 
and t — 2.4 Gyrs. The AWDM galaxy does not follow an 
exponential. There is evidence for a second broad peak in 
the star formation rate around 10 Gyrs ago which is related 



to the enhanced gas accretion of the AWDM galaxy at early 
times. 

The stellar mass assembly history of the 50 3 galaxies 
is shown in Fig. 16. As for the 34 3 simulations, the early 
star formation rate of the central galaxy in the AWDM cos- 
mology is significantly lower than in the ACDM cosmology. 
Due to this difference in the star formation histories, the 
ACDM galaxy is slightly older (half of its stars are older 
than 11.9 Gyrs) than the AWDM galaxy (half of the stars 
are older than 11.2 Gyrs). 

3.3 Convergence Tests: Anticipated Effect of 
Extrapolation to Higher Resolution 

In this section we investigate the effect of numerical resolu- 
tion on the global properties of the galaxies formed in the 
AWDM and ACDM cosmologies. We estimate the degree to 
which convergence has been reached by comparing our 34 3 
and 50 3 runs, which have identical large scale perturbations. 

Table 3 lists some global properties of the 34 3 and 50 3 
galaxies. The only value that increases significantly with res- 
olution for both cosmologies is the radius at which the cir- 
cular velocities reach their peak values. The corresponding 
circular velocity profiles for the baryonic and dark matter are 
shown in Fig. 14. The baryonic matter is more concentrated 
at low resolution. There are a number of possible explana- 
tions for this behaviour. If the gas in the inner parts of the 
galaxy is poorly resolved, star formation leads to a rapid de- 
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Figure 13. Star formation rate of the ACDM and the AWDM 
galaxy versus lookback time. The evolution of the ACDM 
galaxy is close to exponential with a star formation rate oc 
exp (— (t — to)/r) with to = 13.5 Gyrs and r = 2.4 Gyrs. The 
AWDM galaxy shows a strong peak approx. 13.5 Gyrs ago and 
a broad second peak a little more than 10 Gyrs ago. 
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Figure 14. The circular velocity profiles of the 34 3 (upper panel) 
and 50 3 (lower panel) AWDM and ACDM galaxies. The total 
velocity curves are shown by the solid line. The baryonic (dot- 
dashed) and dark matter (dotted) contributions are plotted sep- 
arately. The dynamics of the inner regions (r < 10 kpc) is domi- 
nated by baryonic material, which is more concentrated at lower 
resolution, where as the outer regions of the galaxies are domi- 
nated by dark matter. 
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0.08 


0.53 


0.11 


(1 + z* /2 ) 


3.40 


3.62 


3.07 


3.08 


log\o(Age* (Rmax)/yrs) 


9.971 


10.01 


9.987 


10.01 



Table 3. The physical parameters of the stellar system formed 
in our ACDM and AWDM simulations at the 34 3 and 50 3 reso- 
lution. R^ax is the characteristic radius as calculated from the 
peak in the circular velocity profile, A* is the dimensionless spin 
parameter for the stars, T ro t/T rand is the ratio of the rotational 
kinetic energy to the energy in random motions, M* is the total 
mass within R^, aa; , p*(Rmax) is the density of the stars within 
R^ aI , Zj/2 is the redshift at which half of the stars have formed 
and Age(R^ nax ) is the average age of the stars within H^ nax . 

pletion of gas resulting in artificial pressure gradients which 
drive more gas to the center. There is also the possibility 
that a mismatch of gravitational softening (Plummer) and 
SPH softening (cubic spline) supports unstable regions to 
collapse artificially (see Bate & Burkert 1997). 

The cumulative mass profile is shown in Fig. 15. The 
galaxies produced in the 50 3 simulations extend much fur- 



ther than the lower resolution galaxies due the density cri- 
terion that we use to define the size of the galaxy. 

Fig. 16 compares the stellar mass within 7 kpc, corre- 
sponding to ~ 2r c ff, as a function of lookback time for the 
34 3 and 50 3 simulations in both cosmologies. We see that the 
assembly history is resolution dependent for redshifts z < 2. 
The trend goes in opposite directions for the two cosmolo- 
gies. The AWDM galaxy accretes more mass whereas the 
ACDM galaxy accretes less. The final masses seem to con- 
verge at higher resolution. However, we note that our sim- 
ulations do not include a photo-ionising background which 
would improve the numerical convergence in the case of the 
ACDM cosmology. 

The collapse of objects smaller than the resolution limit 
of the simulations is suppressed at early epochs. Thus, when 
we increase the resolution in the ACDM simulations star 
formation occurs in objects which did not form at the 
lower resolution. However, in the power spectrum of the 
warm dark matter cosmology scales, A < 2.5 Mpc , are in- 
trinsically suppressed. Consequently, we see very little res- 
olution dependence in the star formation history of the 
AWDM simulations as the resolution of the 34 3 simulations 
is A = 0.35 Mpc and therefore well below the cut off in the 
power spectrum. If we were to increase our resolution still 
further we would expect that the early star formation in 
the ACDM simulation to increase whilst very little change 
would be expected in the history of the AWDM simulations. 

If we now compare the effect of resolution on the redshift 
where half of the stars are formed, z*/ 2 , in the two cosmolo- 
gies, we see that it increases with resolution in the ACDM 
cosmology but remains constant for the AWDM cosmology 
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Figure 15. The cumulative mass profiles for the 34 3 (upper 
panel) and 50 3 (lower panel) A WDM (thick) and ACDM (thin) 
simulations plotted for the dark (dot dashed) and baryonic (solid) 
matter separately. Within at least 5 effective radii, the galaxies 
are dominated by dark matter. 
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Figure 16. The evolution of stellar mass within 7kpc (fs 2r c fj) for 
the 34 3 (thin) and 50 3 (thick) simulations in the AWDM (dashed) 
and ACDM (solid) cosmologies. The AWDM simulation appears 
to have converged by the resolution achieved in the 50 3 simula- 
tion, but this is not true of the ACDM simulation. 



(Table 3). Thus if we extrapolate to higher resolution we 
conclude that z* , 2 in the AWDM cosmology will not change 
but the discrepancy between the formation redshifts in the 
two cosmologies would increase and could possibly produce 
a significant observable difference. 



4 DETAILED COMPARISON WITH 
OBSERVATIONS OF ELLIPTICAL 
GALAXIES 

The analysis of the previous sections has shown that the 
simulated galaxies are hot spheroidal stellar systems. In this 
section we present a detailed comparison of the two galax- 
ies formed in the 50 3 ACDM and AWDM simulations with 
observations of giant elliptical galaxies with respect to ob- 
served characteristic photometric and kinematical proper- 
ties, e.g. surface density profile, isophotal deviation from 
perfect ellipses, velocity dispersion, and major- and minor- 
axis rotation. In addition we analyse the line-of-sight ve- 
locity distribution (LOSVD) of both galaxies. The meth- 
ods used to analyse the simulated galaxies closely follows 
the analysis used by observers and is described in detail in 
Naab & Burkert (2003). Only a brief description will be 
given here. We created an artificial image of the simulated 
galaxies by binning the central 35 kpc into 128 x 128 pixels 
and smooth it with a Gaussian filter of standard deviation 
1.5 pixels. Using this image we derived the surface brightness 
distribution and the best fitting Sersic-profile (Sersic 1968) 
E = So exp(— r 1//nsor ), where n sor is the Sersic-index. For 
n sor = 1 the profile is exponential, n scr = 4 parameterises 
the de Vaucouleurs r 1 / 4 profile. The isophotes were analysed 
with respect to the fourth-order cosine deviation, 04, from 
perfect ellipses (Bender et al. 1988). For every projection 
of a simulated galaxy we defined an effective 04-coenicient, 
o4 e ff, as the mean value of 04 between 0.25r e and 1.0r e , 
with r e being the projected spherical half-light radius. For 
the moment we are treating (M/L) as constant for all stel- 
lar particles regardless of when they were created, or the 
metallicity of the gas from which they originated. 

The characteristic ellipticity e e ff for each projection was 
denned as the isophotal ellipticity at 1.5r e . The central 
velocity dispersion <r was determined as the average pro- 
jected velocity dispersion of all stellar particles inside a pro- 
jected galactocentric distance of 0.2r e . The characteristic 
rotational velocity along the major and the minor axis were 
the projected rotational velocities at 1.5r e and 0.5r e , respec- 
tively. 

We analysed the surface brightness distributions of 100 
random projections of the ACDM and AWDM galaxy using 
the Sersic-profile. The ACDM galaxy has a half-mass radius 
of r e g « 3.2 kpc and a mean Sersic-index of n scr = 1.8 ± 0.1 
which is close to the observed lower limit for intermediate- 
mass ellipticals at the same effective radius (Caon et al. 
1993). The AWDM galaxy shows a more exponential profile 
with a mean Seric-index of n sor = 1.5 ± 0.15 at a smaller 
effective radius of r e g ~ 2.2 kpc again consistent with lower 
observational limits (see Fig. 23) . Interestingly, the observed 
trend for more exponential galaxies to have smaller effective 
radii is followed. Combined with the ellipticities and isopho- 
tal deviations derived in the following section the photomet- 
ric properties of both galaxies agree with observed ellipticals 
(see Caon et al. 1993). 

In Fig. 17 we show the global analysis of the two galax- 
ies resulting from the 50 3 ACDM (filled boxes) and AWDM 
(open triangles) simulations seen along 100 random lines- 
of-sight. The observations of real elliptical galaxies are indi- 
cated by plus signs (data kindly provided by Ralf Bender). 
In the upper left plot of Fig. 17 the projections of the ACDM 
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Figure 17. Kincmatical and photometric properties for 100 random projections of the simulated AWDM and ACDM galaxies in the 
50 3 simulation. Open triangles show the values for the AWDM galaxy while filled boxes indicate the values for the ACDM galaxy. 
Observed values are indicated by plus signs, arrows indicate upper limits (Data kindly provided by Ralf Bender). The errors were derived 
by statistical bootstrapping. Upper left panel: Effective ellipticity e e ff of the galaxies versus effective isophotal shape a4 e ff. Upper right 
panel: Rotational velocity over central velocity dispersion (f ma j/<ro) versus e c g . The dashed line shows the theoretical curve for an oblate 
isotropic rotator. Lower left panel: Anisotropy parameter (v/o~)* versus a4 e g . Lower right panel: Amount of minor axis rotation versus 
a4 c ffwith i> ma j an v m i n being the maximum velocity along the major and minor axis, respectively. 



galaxy in the a4 c ff-e e ff plane are shown. The galaxy looks al- 
most round with a maximum ellipticity of 0.2 and a peak in 
the ellipticity distribution around e = 0.1. It has predomi- 
nantly disky (a4 c ff > 0) deviations from elliptical isophotes 
and there is a trend for more disky projections to have higher 
ellipticities. The AWDM galaxy is more elongated and the 
distribution of ellipticities peaks at e = 0.2 with a maxi- 
mum of e = 0.35. The isophotes are either slightly boxy or 
disky. Taking the errors of the 04 analysis into account the 
isophotal shape scatters around zero. The area covered by 
the projections of both galaxies in the a4 e fr-e e ff plane is con- 
sistent with the observed distribution. However, the simula- 
tions do not reproduce galaxies with e e s > and a4 e s > 2 
or aAes < —I. 

With respect to rotation, the ACDM galaxy is a slow 
rotator with (v ma j/co) < 0.2 (upper right plot in Fig. 17) 



and is likely to be flattened by anisotropic velocity disper- 
sions as the anisotropy parameter (v m aj / '(To)* is in general 
well below 0.7. As this galaxy shows predominantly disky 
isophotes it occupies a region in the (v ma j / 'co)* -a^es plane 
where no real galaxy is observed (lower left plot in Fig. 17). 
Additionally, in the (iWj/0"o)-e c ff-plane the projections 
cover the area of observed massive boxy ellipticals which 
are much more luminous than the galaxies presented here. 
In contrast, the AWDM galaxy has projected (v ma j/ao)- 
values higher than 0.5, which is consistent with models 
for rotationally supported galaxies (dashed line in upper 
right plot of Fig. 17). The galaxy also is fairly isotropic 
with (vmaj/cro)* scattering around zero and agrees well 
with observed isotropic galaxies with small deviations 
from elliptical isophotes. None of the projections of the 
simulated galaxies falls in the area of boxy anisotropic 
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Figure 18. Analysis of the LOSVD of the 50 3 ACDM galaxy as 
measured along a slit aligned with the major axis of the moment 
of inertia tenor of the main stellar body (open squares) . The fitted 
local velocity i>g t , local velocity dispersion <Tfj t , u ro t / fflt j ^3 > and 
/14 are plotted versus radius. The lines show the smoothed profiles, 
respectively. The effective radius of the galaxy, r e g ~ 3.2 kpc is 
indicated by the arrows. The individual error-bars were derived 
by bootstrapping. 



or disky isoptropic ellipticals. Both simulated galaxies 
show minor-axis rotation (lower right plot of Fig. 17). The 
ACDM galaxy however disagrees with observations as it 
shows disky isophotes at the same time. Observed disky 
ellipticals show only weak minor axis rotation. 

To analyse the LOSVDs of the simulated galaxies in 
more detail we placed a slit with a width of 1.5 kpc and a grid 
spacing of 0.5 kpc along the apparent long axis of each pro- 
jected remnant. Thereafter all particles falling within each 
grid cell were binned in velocity along the line-of-sight. The 
line-of-sight velocity profiles for each bin along the grid were 
then parameterised by a Gaussian plus third- and fourth- 
order Gauss-Hermite basis functions (van der Marel & Franx 
1993; Bender et al. 1994; Bendo & Barnes 2000). Figs. 18 and 
19 show the kinematic parameters (ufH, Vfit, hz,hi) along the 
long axes of the ACDM and AWDM galaxy, respectively. hz 
and hi are the amplitudes of the third- and fourth order 
Gauss-Hermite functions. For /1.3 = and /14 = the result- 
ing velocity profile is a Gaussian. For asymmetric profiles 
with the prograde (leading) wing steeper than the retro- 
grade (trailing) one, hz and vm have opposite signs. When 
ffit and /13 have the same sign, the leading wing is broad 
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Figure 19. Same as Fig. 18 but for the 50 3 AWDM galaxy. The 
effective radius of the galaxy, r e g ~ 2.2 kpc is indicated by the 
arrows. 



and the trailing wing is narrow. LOSVDs with /14 > have 
a 'triangular' or peaked shape, here the distribution's peak 
is narrow with broad wings. Flat-top LOSVDs have hi < 0, 
where the peak is broad and the wings are narrow. 

The ACDM galaxy shows very little rotation inside 
its effective radius (Fig. 18). The velocity dispersion inside 
r e g is flat and falls off slowly at larger radii therefore 
Wflt/cfit never exceeds 0.2. The asymmetry of the LOSVD 
parameterised by hz is anticorrelated with vat as it is 
observed for real ellipticals (Bender et al. 1994). The local 
correlation between /13 and Vfu/cfit inside e e g for 100 ran- 
dom projections of the simulated galaxy is shown in upper 
plot of Fig. 20. hz and «fi t /o"fl t agree well with observed 
data (indicated by the shaded area) for v^t/vm < 0.2 and 
are in general anticorrelated. This is also reflected in a 
negative effective hz value, /i3 c ff, (lower plot in Fig. 20) for 
almost all projections, where /i3 c ff is defined as the mean 
value of hz inside one effective radius (see Bender et al. 
1994). The points follow the observed correlation between 
/i3 e ff and v TO t/a a . 

The AWDM galaxy, as we have already seen in Fig. 17, 
shows significant rotation (Fig. 19) with the rotation veloc- 
ity flattening out at around 100 km/s beyond the effective 
radius. The velocity dispersion profile is more peaked than 
for the ACDM case. This results in local va t /est as high as 
0.6. The parameters hz and /14 are consistent with zero in- 
side r e ff . Analysed for 100 projections the local value of hz 
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Figure 20. Upper plot: Local correlation between I13 and vat/ao 
for the ACDM galaxy. The observed data points cover the region 
indicated by the shaded area. The large error bar indicates the 
errors for v^t/ao > 0.2. Lower plot: Global correlation between 
h3 e ff and va t /oo The observed correlation (Bender et al. 1994) is 
shown by a straight line. The maximum spread around this cor- 
relation is indicated by the shaded area. The dotted lines indicate 
the zero line. 

does not correlate with vat /cat which is in contradiction to 
observations (upper plot of Fig. 21). h3 e s (lower plot in Fig. 
21) tends to stay positive and does not follow the observed 
trend, however it is broadly consistent with observations. 

In summary, the AWDM galaxy can be interpreted as 
a modestly rotating isotropic hot stellar system. Its global 
photometric and kinematical properties are consistent 
with observations of isotropic ellipticals. The galaxy does, 
however, not show any strong sign of the fine-structure 
that would be typical for rotationally supported elliptical 
galaxies. Its isophotal shape does not deviate significantly 
from being elliptical and there is no evidence that the 
LOSVDs show a correlated deviation from a Gaussian. This 
may be due, at least in part, to an insufficient number of 
particles in the simulation. 

The ACDM galaxy shows signatures for being a 
combined disk/bulge system (except for the very low net 
rotation rate). The photometric decomposition (Fig. 5) 
presents evidence for a disk-like sub-component. This 
is supported by the predominantly disky shape of the 
galaxy. The anticorrelation between /13 and v/a provides 



Figure 21. Same as Fig. 20 but for the AWDM galaxy. 

further kinematical evidence for a combined disk-bulge 
system (Bender et al. 1994; Naab & Burkert 2001b). From 
the rotation curves, the round shape and the anisotropic 
kinematics of the galaxy it becomes clear that the "disk" 
can contribute only a small fraction to the total mass of 
the system. However, the ACDM galaxy rotates too slowly 
to be consistent with observations and there is no sign of 
a disk in the kinematic decomposition. It forms a disky 
anisotropic system which is not observed. 

Due to the uncertain mass-to-light ratios of the 
simulated galaxies we use the virial theorem to define a dy- 
namical masses M dyn = r e scro/G using the half- mass radius 
r e ff and the central velocity dispersion o~o of every projected 
galaxy. This mass is used to define a dynamical surface den- 
sity, Edyn = Afd yn /r;?ff, to create a simulated Fundamental 
Plane where the effective surface brightness SB e is replaced 
by Edyn- In the same way we derive dynamical masses 
for observed ellipticals using the observed effective radius 
and velocity dispersion. Fig. 22 shows a comparison of 
observed ellipticals with the simulated ACDM and AWDM 
galaxies. Both galaxies are consistent with observations. 
The AWDM galaxy has a smaller dynamical mass and has 
a higher effective surface brightness than the ACDM galaxy. 

Figure 23 shows a comparison of the mean global pro- 
jected properties of the simulated galaxies to the global 
properties of giant elliptical galaxies. The maximum spread 
in the data and in the observations is indicated by the er- 



16 Lisa J. Wright, Jeremiah P. Ostriker, Thorsten Naab & George Efstathiou 



1Q' 



4 



1D 1n ro 



|0 9 
10 



I 1 1 1 1 I 




■ ■ ■ ■ ■ ■ 



a. 



I ' ' ' I 



o 



O 9 O 



c 8 



O CP o 




100 



Figure 22. Dynamical mass M(j yn = r e sai/G, effective radius 
r e ff , and dynamical surface- brightness, E dyn , versus central veloc- 
ity dispersion <jq. Black dots represent projections of the ACDM 
, grey dots of the A WDM galaxy, respectively. For the simulated 
galaxies r e g is the projected spherical half-mass radius. 



ror bars. It becomes clear from this figure that we produced 
"ordinary" ellipticals in our simulations which compare very 
well to the average properties of giant ellipticals. 

The low number of stellar particles within the half-mass 
radii of the simulated galaxies (« 7000) makes it too early 
to draw any final conclusion on the real structure of inter- 
mediate mass ellipticals in AWDM or ACDM cosmologies. 
Especially as the ACDM simulations are not yet numeri- 
cally converged. Relaxation effects in particular might signif- 
icantly influence the dynamics of the systems (see Diemand 
ct al. 2003) 

The present resolution is already sufficiently high to 
perform a kinematical and photometric analysis within 
reasonable error limits (all errors have been determined by 
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Figure 23. Mean properties of the simulated ACDM (squares), 
AWDM (triangles) and observed (diamonds) elliptical galaxies. 
The maximum spread in the data is indicated by the error bars. 
For the two simulated galaxies the spread is due to projection 
effects. Observational data for the effective radius (r e g), elliptic- 
ity (e), effective isophotal shape (o4 e ff), rotation velocity (u m aj) 
and velocity dispersion (co) have been kindly provided by Ralf 
Bender. Data for the Seric index n scr have been taken from Caon 
ct al. (1993) and for the LOSVD asymmetries (/i3 e ff) from Bender 
et al. (1994). 



applying the statistical bootstrapping method (see Heyl 
et al. 1994) both for global and local measurements (see 
e.g. Figs. 17, 18). Therefore we are able to trace the real 
structure of the simulated galaxies at the present resolution. 
Increasing the resolution, however, will very likely change 
the detailed internal stellar structure of the simulated 
galaxies whereas the differences will presumably be stronger 
for ACDM (as we have shown in the previous section) than 
for AWDM cosmologies. A preliminary analysis of a higher 
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resolution ACDM supports this conclusion. 



5 CONCLUSIONS 

We have presented the results of ten 34 3 particles and 
two 50 3 particles simulations of the formation of individual 
galaxies in a warm and cold dark matter cosmology with a 
cosmological constant. The sample of low resolution simu- 
lations enabled us to compare the global properties of the 
dark matter haloes and galaxies formed in the two different 
cosmologies. The two galaxies selected for re-simulation re- 
side in a low density environment. Consequently the models 
presented here are are likely to trace the formation of ordi- 
nary intermediate mass giant elliptical or SO/Sa galaxies in 
the field, for example, the Sombrero galaxy (M104). The two 
high resolution simulations have been used to investigate the 
assembly histories of the final galaxies and to compare their 
internal properties with those of real early type galaxies. 

As expected from the suppression of small-scale struc- 
ture in the initial conditions, the AWDM cosmology pro- 
duces fewer low mass dark matter haloes (M^aio < 1O 1O M0) 
at the present epoch compared with the ACDM cosmology 
(Bode et al. 2001). This feature is also reflected in the assem- 
bly histories of the two 50 3 simulations. At almost all red- 
shifts, the ACDM galaxy experiences more minor mergers 
with higher mass-ratios (up to 10:1) than the corresponding 
AWDM galaxy. Low accretion rates, and the absence of any 
major merger event, are expected since the initial conditions 
were selected from low density environments. In addition to 
differences in the frequency of mergers, the internal compo- 
sition of the merging satellites differs in the two cosmologies. 
The AWDM satellites are always more gas rich. At redshifts 
z < 0.3 the star-to-gas ratio is about an order of magnitude 
higher in the ACDM than in the AWDM cosmology. As a 
result, below a redshift of two (when both galaxies have al- 
ready assembled ~ 50% of their final stellar mass) about 
30% of all present day stars in the ACDM galaxy have been 
accreted by mergers and « 20% of the stars have formed 
inside the galaxy. Over the same period the AWDM has ac- 
creted only « 10% of its present day stars and « 40% have 
formed within the galaxy. 

A further difference between the two cosmologies is ev- 
ident in the angular momentum evolution of the galaxies. 
The halo of the ACDM galaxy has higher angular momen- 
tum than the AWDM halo but they both show a similar 
temporal evolution. In contrast, the average angular mo- 
mentum of the stars of the ACDM galaxies does not change 
significantly after a redshift of z — 4 whereas the stars of 
the AWDM galaxy gain further angular momentum until 
z — 1, with little evolution thereafter. Between a redshift of 
z = 4 and z = 1 the AWDM galaxy is more gas rich and 
forms more stars within the galaxy than its ACDM counter- 
part, resulting in a second peak in the star formation rate 
at z ?s 1.5. 

We conclude that the removal of small scale power re- 
duces the frequency of mergers. The reduction in the number 
of massive sub-halo mergers in the 50 3 warm dark matter 
simulation produces a galaxy with significantly higher angu- 
lar momentum at z=0. The increase in the specific angular 
momentum of objects in the warm dark matter simulations 



suppresses the collapse of gas within the simulations, which 
in turn reduces the star formation rate of the central galaxy 
at early epochs. 

We performed a photometric and kinematical decom- 
position of the main stellar systems in the 50 3 simulations. 
Both galaxies are dominated by a hot spheroidal component. 
The global projected properties of the two 50 3 ACDM and 
AWDM galaxies resemble those of real elliptical galaxies. 
The total masses, density profiles, effective radii, ellipticities, 
global values of isophotal shapes and LOSVD asymmetries 
are consistent with the average global properties of giant el- 
liptical galaxies of intermediate masses. This is in contrast 
to properties of very massive, anisotropic, and boxy ellipti- 
cals with large effective radii, typically found in clusters of 
galaxies. As the physical properties of massive ellipticals dif- 
fer from those of ordinary intermediate mass ellipticals with 
respect to inner density-profiles, sizes, isophotal shapes and 
X-ray and radio properties (Bender et al. 1989; Caon et al. 
1993; Faber et al. 1997), we can conclude that they must 
have formed in a denser environment and, therefore had dif- 
ferent formation histories. 

Investigating the AWDM galaxy in detail revealed that 
it appears to be an isotropic fast rotator with only weak fine 
structure both in its isophotal shape (the isophotes are el- 
liptical) and in its LOSVDs (which on average have a Gaus- 
sian shape). The ratio of T rot /T ranc j is a factor of ten higher 
for the AWDM than for the ACDM system. Observed rota- 
tionally supported ellipticals with similar rotational support 
do, however, show on average stronger asymmetries in their 
LOSVDs. The detailed photometric and global kinematical 
properties do agree very well with observations. The half- 
mass radius of stellar population is 2.3 kpc. In contrast, 
the ACDM galaxy shows only weak rotation and appears 
to have anisotropic velocity dispersions. The isophotes are 
predominantly disky and the LOSVDs show the observed 
trend which indicates the presence of a weak disk compo- 
nent embedded in the spheroidal body of the galaxy (Rix 
& White 1990; Bender et al. 1994; Naab & Burkert 2001b). 
However, anisotropic disky systems, like the ACDM ellipti- 
cal, are in general not observed. The half-mass radius of the 
ACDM galaxy is 3.2 kpc. The mean age of the stellar pop- 
ulation of both galaxies is about 10 Gyrs which is in good 
agreement with the ages of early type galaxies. 

Meza et al. (2003) have recently simulated the forma- 
tion of an individual elliptical galaxy in a ACDM cosmology 
which is only a little more massive (by a factor of about 
1.5) than our ACDM galaxy. However, their galaxy does 
not match the global properties of observed giant ellipticals 
as its stellar distribution is far too dense. Numerical resolu- 
tion is unlikely to cause the difference as their initial particle 
masses and softening lengths are comparable to the values 
used for the simulations presented here. Their final galaxy, 
however, has significantly more stellar particles (w 65000) 
than our galaxies (~ 13000). This is due to differences in the 
star formation algorithms. Meza et al. (2003) split their gas 
particles into several stars which can have different masses. 
As a result their mean stellar particle mass is a factor of four 
smaller than our fixed stellar mass. It is not clear in how far 
this difference influences the results, e.g. due to mass segre- 
gation effects in the stellar distribution. A further difference 
is that Meza et al. (2003) make their stars at a rate pro- 
portional to the gas density with a relatively low efficiency. 



18 Lisa J. Wright, Jeremiah P. Ostriker, Thorsten Naab & George Efstathiou 



In the simulations presented here a gas particle is turned 
into a star particle as soon as its density is above a certain 
threshold for more than a dynamical time. As a result the 
gas in the Meza et al. (2003) simulation can collapse to much 
higher densities before stars are formed. 

In addition we omitted feedback processes in our simu- 
lations. Meza et al. (2003) implemented thermal and kine- 
matical feedback. Although one might naively expect that 
the inclusion of feedback would result in less compact ob- 
jects, this goes in the opposite direction to explain the dif- 
ferences between their results and ours. Another possible 
explanation for the differences is the merger history. Our 
galaxies experience only minor mergers with mass ratios up 
to 10:1. The Meza et al. (2003) galaxy undergoes a late ma- 
jor merger with a mass-ratio of 3:1 which could effectively 
drive gas into the centre and convert it into stars. However, 
the star formation history is not significantly different from 
our ACDM galaxy which undergoes no major merger. The 
merger is therefore more likely to influence the dynamics of 
the system rather than have a major effect on the concen- 
tration of the bulk of the stellar population. 

Based on the simulations at two different resolutions in 
a ACDM and AWDM cosmology it becomes clear that this 
investigation is not yet definitive. As the simulations are not 
yet fully numerically resolved (especially the ACDM simu- 
lation) we cannot say with confidence whether ACDM or 
AWDM produces galaxies that most closely resemble real 
elliptical galaxies. In conclusion, we can state that with the 
simple physics included in our simulations, it is possible to 
produce ellipticals in both cosmologies with global proper- 
ties that are in good agreement with observations of interme- 
diate mass giant ellipticals or SOs. However, the combination 
of the detailed properties of our simulated galaxies (which 
are very likely to be resolution dependent), like the shape 
of the LOSVD or the isophotal shape, differ slightly from 
observations of real ellipticals for one combination or the 
other. Future simulations at higher resolution will hopefully 
enable us to determine which, if either, of the two cosmolo- 
gies produce galaxies which match the observations. 
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